Coil-Globule transition of a single short polymer chain - an exact enumeration study 
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We present an exact enumeration study of short SAWs in two as well as three dimensions that addresses the 
question, 'what is the shortest walk for which the existence of all the three phases - coil, globule and the theta 
- could be demonstrated' . Even though we could easily demonstrate the coil and the globule phases from Free 
Energy considerations, we could demonstrate the existence of a theta phase only by using a scaling form for the 
distribution of gyration radius. That even such short walks have a scaling behavior is an unexpected result of 
this work. 
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Proteins are biologically active, long chain heteropolymers 
that exist either in the collapsed (native folded or molten glob- 
ule) phase or in the extended coil phase, depending on their 
environmental conditions (]]]. Simple lattice models iQ] have 
been devised for understanding the role of energy-entropy bal- 
ance in realizing these phases. In fact, even short lattice chains 
consisting of only sixteen monomers or less have been shown 
J^t] to be long enough to throw light on the various phases 
of proteins. These chains are basically Self Avoiding Walks 
(SAW) on a regular lattice, modeling homopolypeptides with 
appropriately parametrized intrachain interactions. 

Long chain molecules, modeled as SAWs, are known to ex- 
ist in one of the three conformational phases fl4|] - namely, ex- 
tended coil phase, theta phase and the collapsed phase - char- 
acterized by distinct values of the universal size exponent 0]. 
For example, in two dimensions, the size exponent v = 3/4 
(extended phase), 4/7 {theta phase) and 1 /2 (collapsed phase) 
respectively. It is again the change in the energy-entropy bal- 
ance that brings about transition from one phase to the other. 

In analogy to the lattice protein models, we may ask what is 
the shortest SAW for which the existence of all these phases 
can be demonstrated, eventhough a clear-cut distinction be- 
tween these phases is ascribed only to asymptotically long 
(N — > oo) walks. Moreover, the possibility of experimen- 
tally observing conformations of a single chain, which is of 
finite length, and locating the coil-globule transition |0] begs 
the question whether the theta phase of a very short chain can 
be identified. 

In this brief report, we first demonstrate coil-globule tran- 
sition for a short SAW consisting of twenty steps or less by 
exact enumeration of its conformations on two dimensional 
as well as three dimensional lattices. Then, using the scaling 
form for the distribution of gyration radius ($], we show that 
theta phase can be clearly identified only for walks having a 
minimum of sixteen steps. The applicability of scaling ideas 
to such short walks is indeed surprising. 

In order to study the thermodynamic properties of SAW, we 
need to assign energy to any given conformation. A standard 
way of doing this for a lattice SAW is to assign a quantum of 
energy, e, to each pair of non-bonded nearest neighbor sites 



that belong to the walk |5|]. So, a conformation with n such 
pairs of sites, also called contacts, has energy E = ne. At any 
given inverse temperature, 0, the canonical Partition Function 
for iV-step SAWs can be written as 
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where nx (N) is the maximum number of contacts an iV-step 
walk can have, and cjv (n) is the number of iV-step SAW con- 
formations with n contacts. It is clear that conformations with 
a given number of contacts will acquire more (less) weights if 
e is negative (positive). Since we are interested in coil-globule 
transition, we can set e = —1 without loss of generality. 
We may rewrite the above Partition Function in the form 



n=nx (-^V) 



Cn{(3) = 



-/3F(n,0) 



n=0 



-f3F{n,(3) = (3n + ln[c N {n)] 



(2) 



(3) 



where F(n, (3) is the Free Energy associated with walks hav- 
ing n contacts. Since cx{n) has a bell-shaped form, for a 
given value of (3, the Free Energy F(n, (3) is a minimum at 
a certain value of n, say np(f3), which could in general be a 
non-integer. 

Another way of estimating np{0) is to write the Free En- 
ergy as an implicit function of the squared end-to-end dis- 



- f3F{n p (r 2 ),l3) = (3Mr 2 ) + ln[c N (n p (r 2 ))} (4) 

where np{r 2 ) is the canonical average of the number of con- 
tacts made by all the conformations whose end-to-end dis- 
tance squared is r 2 , and Cjy(fif}(r 2 )) is the corresponding in- 
terpolated value. The estimate of n^ p at which the Free En- 
ergy is minimum agrees very well with that of uf{(3) defined 
earlier. 

The temperature dependence of fip t p per step, presented 
in Fig[T|for SAWs on Square as well as Triangular lattices, 
clearly suggests a transition from a low temperature phase to 
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FIG. 1: Exact enumeration data np t F /N for two dimensional SAWs 
as a function of (3. (i) Lower one: Square lattice walks of length N = 
16(filled star), 18(inverted open triangle) and 22(shaded filled circle), 
(ii) Upper one: Triangular lattice walks of length N = 14, whose 
ordinates are marked on the right. INSET: n/3,F /N as a function of 
< r 2 >jv (/3) for square lattice walks of length N = 18. 
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FIG. 2: Exact enumeration data np^/N for three dimensional 
SAWs as a function of /3. Open stars correspond to walks of length 
N — 14 on a Cubic lattice (Data shifted up by 0.25). The rest of the 
data correspond to Diamond lattice walks of length N = 16(open 
triangle), 18(open circle) and 20(shaded filled circle) respectively. 



a high temperature phase. In fact, the specific heat associ- 
ated with this transition (C sp .ht = —(Pdnp^/dP) is seen to 
peakat/3pG ~ 0.85, 0.52, 1.85 and 1.35 on square, triangular, 
diamond and cubic lattices respectively. What these phases 
are becomes clear from the inset of this figure which shows 
a monotonic decrease of fip,F P er step as a function of the 
canonically averaged end-to-end distance squared. In other 
words, the low temperature phase corresponds to compact 
walks with large number of contacts whereas the high temper- 
ature phase corresponds to extended walks with small number 
of contacts. We observe a similar transition from low tem- 
perature compact phase to high temperature extended phase 
for three dimensional (Cubic and Diamond lattice) SAWs also 
(see, Figf2]). 

Thus, fij3,F per step , or equivalently the canonically av- 
eraged end-to-end distance squared, seems to be a physical 
quantity that is responsive to the coil-globule transition. The 
latter has the advantage of being proportional to the measur- 
able average radius of gyration. But, they provide no evidence 
for the existence of theta phase. This leads to the question of 
how the geometric parameters characterizing a SAW - namely, 
number of steps in the walk (AT), number of contacts made (n) 
and the radius of gyration or the end-to-end distance(r 5 or r)- 
can be put together so as to check whether theta phase for very 
short walks is identifiable. 

Lhuillier |01 first proposed a scaling form for the distribu- 
tion of gyration radius and showed that it is related to the Flory 
free energy for SAW. In fact, radius of gyration is a better ge- 
ometrical descriptor of a SAW than its end-to-end distance. 
For example, in the globule or coil regime, smaller (larger) ra- 
dius of gyration implies larger (smaller) number of contacts; a 



similar correlation between end-to-end distance and the num- 
ber of contacts is not so obvious because conformations with 
different end-to-end distances could all have the same radius 
of gyration. In the stretched regime, however, radius of gy- 
ration and end-to-end distance provide equivalent description 
because of very few contacts. 

Subsequent studies [8] along these lines suggest that the 
most relevant order-parameter, t, for studying the various 
phases of SAW is a certain power of the monomer density, 
namely, t = p 1 /^ -1 ) where p = N/r^ and d is the dimen- 
sionality of space. The distribution of t for an iV-step SAW 
has a scaling form (H] given by 

P N (t) ~ t c exp{-[ ai t + a 2 t 2 + a 3 t- q }); q = V ~ X ' d (5) 

l — i/ 

where c is a universal exponent, related to the susceptibil- 
ity exponent 7, having values —37/32 and ~ —1.13 in two 
and three dimensions respectively. The exponent q = 1 and 
~ 0.617 in two and three dimensions respectively. The pa- 
rameters, a 1 ,a 2 and (13, are /3-dependent and, in particular, 
the value of (3 at which the parameter ai becomes zero corre- 
sponds to the inverse 9 temperature, say fig. 

We have shown in Fig [3] the distributions, Pjv(t), corre- 
sponding to inverse temperatures = 0.0, 0.4, 0.8 and 1.5 re- 
spectively, obtained from exact enumeration data for 23-step 
walks on a Square lattice. By fitting the above scaling form, 
Eq.©, to the histogram data obtained for various values of 
j3, we find that the parameter a\ = for j3 = f3g ~ 0.33. 
Repeating this procedure, we obtained (3g values for N = 14 
to 23 and have presented them in FigJH they stabilize around 
an average value ~ 0.32 for iV > 16. Similarly, the (3g val- 
ues obtained for walks of length in the range 13 to 22 on a 
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FIG. 3: Distributions of order parameter t = (N /r 2 .) 2 for Square 
lattice walks of length N — 23 corresponding to /3 — 0.0, 0.4, 0.8 
and 1.5 (respectively from left to right). Solid lines are the scaling 
form Eq. ([5} fit to the the data. 
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FIG. 4: /3$ values at which the parameter ai = (Eq.lO) for various 
values of N in the range 14 to 23 for Square lattice (lower figure), 
and in the range 13 to 22 for Diamond lattice (upper figure). 

Diamond lattice (the upper curve in Fig|4]l are also seen to 
stabilize around an average value ~ 0.65 for N > 15. This 
implies that a theta phase could be identified for even for such 
short SAWs. 

It may be noted that the fit (solid lines in FigO is not as 
good in the large t or 'globule' regime as it is in the small t or 
'coil' regime. We found this to be the case for Diamond lat- 
tice walks also. The reason could be that the walk is too short 
for the scaling form given by Eq.((5]l to hold good in the entire 



range. Moreover, the average values of (3g, obtainable from 
the data presented in Fig|4j are different from the correspond- 
ing values obtained for asymptotically long walks (N — > oo). 
Nevertheless, they provide clear evidence for the existence of 
a theta phase where the parameter a% vanishes. 

In summary, using exact enumeration data, we have demon- 
strated the existence of theta phase besides the usual coil and 
globule phases even for short SAWs (N < 23) in two as 
well as three dimensions. This is analogous to realizing the 
various phases of proteins by exactly enumerating short lat- 
tice (HP) chains with appropriately parametrized interactions. 
Even though the coil and the globule phases could easily be 
demonstrated from Free Energy considerations, the existence 
of a theta phase could be demonstrated only by using a scaling 
form for the distribution of the gyration radius. That scaling 
ideas work for such short walks is a surprising result. 
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